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ABSTRACT 

Luminous extra-galactic water masers are known to be associated with active galactic nuclei 
and have provided accurate estimates for the mass of the central supermassive black hole and 
the size and structure of the circumnuclear accretion disk in nearby galaxies. To find water 
maser systems at much higher redshifts, we have begun a survey of known gravitationally 
lensed quasars and star-forming galaxies. In this paper, we present a search for 22 GHz (rest 
frame) water masers toward five dusty, gravitationally lensed quasars and star-forming galax- 
ies at redshifts between 2.3 and 2.9 with the Effelsberg radio telescope and the Expanded Very 
Large Array. Our observations do not find any new definite examples of high redshift water 
maser galaxies, suggesting that large reservoirs of dust and gas are not a sufficient condition 
for powerful water maser emission. However, we do find the tentative detection of a water 
maser system in the active galaxy IRAS 10214+4724 at redshift 2.285. Our survey has now 
doubled the number of gravitationally lensed galaxies and quasars that have been searched 
for high redshift water maser emission. We also present an updated analysis of the high red- 
shift water maser luminosity function that is based on the results presented here and from the 
only cosmologically distant (z > 1) water maser galaxy found thus far, MG J0414+0534 at 
redshift 2.64. By comparing with the water maser luminosity function locally and at moder- 
ate redshifts, we find that there must be some evolution in the luminosity function of water 
maser galaxies at high redshifts. By assuming a moderate evolution [(1 + z) 4 ] in the wa- 
ter maser luminosity function, we find that blind surveys for water maser galaxies are only 
worthwhile with extremely high sensitivity like that of the planned Square Kilometre Array 
(Phase 2), which is scheduled to be completed by 2020. However, instruments like the EVLA 
and MeerKAT will be capable of detecting water maser systems similar to the one found from 
MG J0414+0534 through dedicated pointed observations, providing suitable high-redshift tar- 
gets can be selected. 

Key words: gravitational lensing: strong - masers - radio lines: galaxies - galaxies: nuclei: 
starburst 



1 INTRODUCTION 

The 6i6— 523 transition of water (emitted rest-frequency of 
22.23508 GHz) shows spectacularly luminous extra-galactic maser 
emission, with iso tropic lin e luminosities of up to 23 000 Lq (for 
a review see, e.g., lLol l2005h . The maser action results from colli- 
sional excitation of water molecules in extremely dense gas clouds 
and requires particle densities of 10 7 to 10 11 cm -3 and gas tem- 
peratures > 300 K. All interferometric studies of powerful water 



masers (> 10 Lq) found thus far have always shown them to be as- 
sociated with active galactic nuclei (AGN). High-resolution studies 
with very long baseline interferometry (VLBI) of ~ 10 of these sys- 
tems have shown that the water masers are associated with either 
the circumnuclear accretion disk within a few parsecs of the cen- 



tral supermassive black hole jMivoshietai1ll995l :l Herr nstein et alj 
1 19991 : kondratko et alj [20081 : iReid et al.1 120091) or within the in- 
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ner part (< 30 p c) of the relativistic radio jets that a re ejected 
from some AGN dClaussen et al.lll998l : |Peck et alj|2003h . The high 
surface brightness and proximity to the central engine make wa- 
ter maser systems a powerful tool for studying the properties of 
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AGN. For example, water masers have been used to measure the 
mass of the central supermassive black hole i Mivoshi et al.ll 19951 ; 
iGreenhill et al]|2003al ; iReid et al.ll2009l ; iKuo et al.l[2oioir the size 



and structure of the circumnuclear accretion disk ( Greenhill et al 



2003al ; lArgon et alJl2007l ; lHumphre"vse t al. 2008; Kondratk oet alJ 
2008), the shock speeds and densities of radio jets jPeck et al.l 
2003J), a nd to determine accurate geometrical dist ances to the host 



galaxies jHerrnstein et alJI 19991 ; IBraatz et al.| [2010) - all of which 
are competitive with the results obtained using other methods. 

Surveys for water maser systems at low redshifts (z ^ 0.06) 
have tended to focus on Seyfert 2 an d LINER galaxies (e.g. 
IGreenhill et ai]|2003bl ;l Braat z et al. IHooi). This is because the ex- 
pected edge-on orientation of the circumnuclear accretion disk for 
these galaxies should result in a larger path length of velocity- 
coherent molecular gas which is needed for the maser ampli- 
fication. The success rate of these surveys is of order ~5 per 
cent, with a round ~130 water maser galaxies currently known 
locally (e.g. iBraatz et all 2008; and unpublished). Surveys have 
also targeted nearby spiral an d starburst galaxies with bright far- 
infrared (FIR) flux-den isties faenkel etall 120051 ; ICastangia et al.l 
l2008l;ISurcis et al.l2009h . These surveys have found 15 water maser 
galax ies and have a detection rate of 23 per cent jSurcis et al.l 
l2009h . A snapshot survey of 611 luminous active and non-active 
galaxies in the loca l Unive rse (v ^ 5000 km) was carried out by 
IBraatz & GugliuccH ( |2008l) . They found only 8 water maser sys- 
tems (a detection rate of about ~ 1 per cent), which demonstrates 
the rarity of these objects. Finding water maser galaxies at higher 
redshifts has the potential to constrain the nature of dark energy 
through the accurate measurement of geometrical distances. How- 
ever, due to the limited sensitivities and the small range of frequen- 
cies covered by current radio telescopes, surveys fo r high redshift 
water maser galaxies had yielded non-detections JBennert et al.l 
2009), with the notable exception being the detection of the most 
lumi nous water maser system kno wn in a type 2 quasar at redshift 
0.66 jBarvainis & Antonuccill2005h . 

To overcome the sensitivity limitations we have started a 
search for water masers from known gravitationally lensed quasars. 
Observing only gravitationally lensed active galaxies has the ad- 
vantage of pre-selecting a population of distant AGN (z > 1) and 
of using the magnification provided by the foreground gravitational 
lensing galaxy to increase the observed flux density of any water 
masers in the background AGN. Such lensing magnifications can 
range up to hundreds, depending on the alignment of the lens and 
the background source. The gravitational lensing technique has al- 
ready been used successfully to study the interstellar medium and 
molecular gas content of high redshift quasars through the d etec- 
tion of, for example, carbon monoxide jBarvainis et al.|[l994h . We 
found a luminous water maser in the first gravitationally lensed 
quasar we observed, in the type 1 quasar MG J0414+0534 at red- 
shift 2.64 ( Impell izzeri et alj2008h . which is by far the most distant 
object to show water vapour emission. The maser line from MG 
J04 14+05 34 is broad with a full width at half maximum (FWHM) 
of ~ 45 km s _1 , is blueshifted by —300 km s _1 from the systemic 
velocity of the quasar and has an (unlensed) apparent isotropic lu- 
minosity of ~8300 1/0. Our initial hypothesis, based on spectra 
obtained with the Expanded Very Large Array (EVLA) and the Ef- 
felsberg radio telescope, was that the maser originates in the rela- 
tivistic jet of this quasar as it interacts with a molecular cloud lying 
close to the supermassive black hole. Follow-up VLBI and moni- 
toring observations are currently being analyzed to investigate this 
hypothesis and will be presented in forthcoming papers. 

We have now begun a survey of other distant gravitation- 



ally lensed quasars, and have also extended our search to lensed 
starbursting galaxies, using the Effelsberg radio telescope and the 
EVLA. The aim of this survey is to find additional detections of 
high redshift water masers, the results of which are presented in 
this paper. In the case of the only currently confirmed high redshift 
water maser galaxy, MG J0414+0534, we know that the source is 
highly dust reddened and hosts a powerful AGN that is emitting a 
relativistic jet. This suggests that luminous water masers that are 
associated with AGN are more likely to be found in high redshift 
galaxies with large amounts of dust and/or powerful radio jets. To 
investigate the first of these possible effects, we have limited the 
sample studied in this paper to include only radio-quiet quasars and 
star-forming galaxies with previously measured strong far-infrared 
(FIR) emission or molecular CO line emission. 

For all calculations we adopt an Qm = 0.3, 
SIa = 0.7 spatially flat Universe, with a Hubble constant 
of -ffo = 70 kms -1 Mpc -1 , and a solar luminosity of 
L = 3.939 x 10 26 W. 



2 TARGETS 

We targeted and prioritized those gravitational lens systems with 

i) background AGN or star-forming galaxies at redshifts between 
2.29 and 2.87 (due to the bandpass of the receiver that was used), 

ii) had high lensing magnifications and iii) showed previous detec- 
tions of molecular emission (e.g. CO, HCN, etc.) or a large FIR 
luminosity. We also limited our sample to those objects that were 
observable with the Effelsberg radio telescope and the EVLA. The 
resulting sample is not statistically complete. 



2.1 RXJ0911+0551 

RX J091 1+0551 is co mprised of four images of an X-ray d etected 
quasar at redshift 2.80 jBade et al.lll997l ; lBurud et al.ll998h . which 
is gravitationally l ensed by a group of galaxies at redshift 0.769 
dKneib et alj|2000h . The total magnification of the quasar i s ~ 22 . 
CO (3-2) has been tentatively detected by lHainline et all (2004) 
and the large sub-mm lumin osity of the lensed quasa r suggests a 
cold dust mass of ~10 8 M Q jBarvainis & Ivisonll2002l) . 



2.2 IRAS 10214+4724 

IRAS 10214+4724 is a highly magnified a nd dust obscured gravi- 
tation ally lensed quasar at redshift 2.285 jRowan-Robinson et al.l 
Il99lh . This system has four lensed images that are formed 
by a massive foreground lensi ng galaxy at redshift 0.893 
jLacv. Rawlings & Serjeantl Il998l) . The m agnification of the 
three merging lensed images is < 100 (e.g. iBroadhurst & Lehan 
1995). The lensed source is extremely lu minous in the FIR 
( Rowa n-Robinson et al.l Il99ll) and sub-mm jBarvainis & rvisonl 
l2002h. and optical po larimetry found evidence of scattered light 
(Lawrence et al. Ill993h. A la rge quantity o f molecular gas from 
CO jBrown & Vanden Bou3 1 19911 1 19921 ; ISolomon et al] 1 19921; 



iDownes et aT1ll995l), HCN jVanden Bout e t al. 2004) and neutral 
carbon ( Weifi et all 120051) has also been reported. Finally, near- 
infrared imaging spectroscopy found evidence for both circumnu- 
clear star formation and an AG N from narrow (bu t spatially ex- 
tended) and broad Ha emission jKroker et alJI 19961) . All of these 
data are consistent with the central engine of the AGN being ob- 
scured and there being nuclear driven star-formation in a region 
rich in molecular gas. 
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2.3 H 1413+117 - The Cloverleaf 

H 1413+117 has four lensed images of a radio-quiet q uasar at red - 
shift 2.560 that are in the shape of a cloverleaf dMagain et al.ll988l) . 
The optical emission from the four lensed images is so strong that 
attempts to measure the redshift of the lens have not been suc- 
cessful. The total magnification of th e background quasar is esti- 
mated to be ~ 10 dKavser et alj |T990*). The system has a plethora 
of m olecular emission and absorption det ections; for example, CO 
(3-2; lBarvainis et all 1991 4-3, 5-4 7-6; lBarvainis et alj 19971 and 
5-6, 8-7, 9-8; iBradford et al]|2009h HCN (1-0; ISolomon et all 
l2003h and CN (3-2; iRiechers et all |2007» . The quasar is also 
ultra-luminous in t he sub-mm and has an estimated dust mass 
of ~ 4 x 10 8 Mrs, ( iBarvainis & Ivisonl l2002h . Previous observa- 
tions searching for maser activity failed to detect any water maser 
emission from this quasar down to an unlensed isotropic lumi- 
nosity limit of 46 00 Lq, with a 10 kms -1 spectral resolution 
dWilneretal] 19991) . 



2.4 MS 1512-cB58 

MS 1512— cB58 is a star-forming Lyman break galaxy at redshift 
2.727 which is gravitationally lense d by t he massive cluster MS 
1512+36 at redshift 0.37 free et alj[l996h . The galaxy is lensed 
into a gravitational arc by th e cluster and is es timated to have a high 
total magnification of > 50 dSeitz et alj|l998l) . CO (3-2) has been 
detected from the galaxy dBaker et alj|2004h~ which along with the 
large potential maser magnification led to MS 1512— cB58 being 
included in our sample. 



2.5 SMM J16359+6612 

SMM J16359+6612 is composed of three lensed images of a pair 
of merging galaxies at redshift 2.516. The lens is the massive clus- 
ter A 2218 at redshift .18, and provides a total magnification of 
~45 dKneib et alj2004h. The mergin g lensed galaxies show strong 
sub-mm emission ( Kneibetal. 2004) and a number of CO detec- 
tions (transitions 3-2, 4-3, 5-4, 6-5, 7-6; e.g.. lWeifietal]|2005l ; 
iKneib et alj|2005l) . SMM J16359+6612 was recently observed for 
water maser emiss i on wit h the E VL A by Ed monds et ai] d2009l) and 
Wagg & Momjian (2009). However, no water vapour was detected 
down to an unlensed luminosity of ~ 5400 Lq within a 10 km s _1 
channel. 



3 OBSERVATIONS 

In this section, our observations with the Effelsberg 100 m radio 
telescope and subsequent follow-up observations with the EVLA 
of tentative water maser detections are presented. 



3.1 Effelsberg observations and analysis 

The spectroscopic observations with the Effelsberg 100 m radio 
telescope were carried out with the 5 cm dual-polarization HEMT 
receiver at the primary focus. This receiver is sensitive between 
5.75 and 6.75 GHz, which for the detection of 22.2 GHz water 
masers corresponds to a redshift range of 2.29 to 2.87. The pri- 
mary beam of the telescope at these frequencies has a FWHM of 
~120 arcsec, which is easily large enough to encompass the angu- 
lar extent of the gravitationally lensed images. 

The data were taken using the position-switching mode where 



a 2.5 min on-source scan was immediately followed by a 2.5 min 
off-source scan. This off-source scan was subtracted from the on- 
source scan to remove the contribution of the sky background 
and the instrument from the on-source measurements. Observa- 
tions of standard flux-density ca librators were also taken to deter- 
mine the antenna gain (3C295; lott et alj|l994l) . The spectra were 
formed usin g the 16 384-chann el Fast Fourier Transform Spectrom- 
eter (FFTS ; (Klein et all 120061) . Both the 20 and 100 MHz band- 
widths that are available with the FFTS were used during these 
observations. A summary of the integration times, observing dates 
and bandwidths used for each lens system is given in TableQ] 

The spectra were analysed using the CLASS77 package within 
GILDAS. Each scan was inspected for any strong gain variations 
or radio frequency interference (RFI), with any spurious scans or 
frequency ranges flagged. The scans were then averaged to form 
a spectrum of the lens system for each observing epoch. A low- 
order polynomial was then fitted to the spectrum to subtract any 
continuum emission from the lensed images. For those lens sys- 
tems that were observed on more than one occasion, an average 
spectrum weighted by the integration times was also produced. We 
smoothed the spectra using Hanning smoothing to a spectral reso- 
lution of 10 kms -1 channel -1 . Where we found a possible water 
maser line, a Gaussian line profile was fitted to the unsmoothed 
spectrum to determine the FWHM, the integrated line flux, the line 
peak flux density and the line velocity. 



3.2 Expanded Very Large Array observations 

Our spectroscopic observations found tentative detections of water 
maser lines from the gravitational lens system IRAS 10214+4724 
(see Section [4] for details). To confirm these detections indepen- 
dently, we carried out interferometric spectroscopic observations 
with the EVLA. These observations would also spatially match any 
water maser emission with the lens system. This is the same strat- 
egy that was used to confirm the d etection of a water mase r from 
the lensed quasar MG J0414+0534 dlmpellizzeri et alj|2008l) , 

IRAS 10214+4724 was observed at a central frequency of 
6.7675 GHz with the EVLA on 2008 May 30 and 31 using 15 
antennas that had recently been upgraded with the new C-band 
receivers. The data were taken through a total bandwidth of 6.25 
MHz and split into 64 channels. This gave a spectral coverage 
of 277 kms" 1 and a spectral resolution of 4.3 kms -1 per chan- 
nel. IRAS 10214+472 4 was known to be a weak radio source 
dLawrence et alj Il993h . therefore the observations were phase- 
referenced by observing the calibrator JVAS J095 8+474 every ~ 15 
minutes to determine the phase and amplitude solutions. Observa- 
tions of another nearby continuum source, JVAS J 103 3+4 12, were 
also taken throughout the run to verify the calibration process. The 
flux-density calibration was determined with 3C 147 and 3C 286 
assuming 6.7 GHz flux-densities of 5.858 and 6.051 Jy, respec- 
tively. IRAS 10214+4724 was observed for ~ 8.5 h in total. The 
data were reduced in the standard way using CASA. A data cube 
was created using all of the channels, except for the first and last 
five channels. The noise per channel was 0.57 mjy beam" 1 and the 
noise in the resulting continuum map was 78 /j,Jy beam -1 . 



4 RESULTS 

We now present a brief description of the radio spectra obtained for 
each gravitational lens system. Isotropic line luminosities, and up- 
per limits in the case of non-detections, have been calculated using 
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Table 1. Observing logs for the Effelsberg radio telescope observations of gravitationally lensed galaxies. 



Lens system 


RA 


Dec 


Redshift 


Date 


Time 


Centre frequency 


Velocity 


range 


Bandwidth 




(J2000) 


(J2000) 






(min) 


(GHz) 


(km s 


- 1 ) 


(MHz) 


RXJ091 1+0551 


09:11:27.500 


05:50:52.00 


2.793 


2008 July 04 


152.5 


5.8617 


-2018 : 


2207 


100 


IRAS 10214+4724 


10:24:34.599 


47:09:11.00 


2.285 


2007 November 12 


209.1 


6.7678 


-267 : 


252 


20 










2008 March 14 


116.1 




-313 : 


304 


20 










2008 March 17 


215.7 




-266: 


252 


20 


H 1413+117 


14:15:46.300 


11:29:43.00 


2.561 


2008 August 15 


283.7 


6.2440 


-1898 : 


2113 


100 


MS 1512-cB58 


15:14:22.219 


36:36:24.79 


2.720 


2008 May 27 


103.1 


5.9772 


-2014 : 


2216 


100 










2008 August 15 


87.8 




-1884: 


2154 


100 


SMM J16359+6612 


16:35:52.501 


66:12:15.01 


2.520 


2008 May 26 


332.5 


6.3167 


-1846 : 


2022 


100 



L H2 o _ 1 0.023 JSdu D\ 
Lq m 1 + z [Jykms -1 ] [Mpc 2 ] ' ' 

where m is the estimated lensing magnification, z is the redshift of 
the background source, J S dv is the integrated line profile, Dl is 
the luminosity distance and Lh 2 o is the isotropic luminosity. We 
quote luminosities after correcting for the estimated magnification 
of the lens. However, these values should be taken with caution 
since the actual magnification of any given maser region can be 
calculated only from high-resolution imaging of the masers (to de- 
termine their position relative to the lensing galaxy) and lens mod- 
elling. This is particularly important for those cases with extended 
lensed emission where there can be a large differential magnifi- 
cation over the extent of the lensed images. All of the luminosity 
upper limits are for rectangular maser lines of width 10 km s - 1 . 

4.1 RX J0911+0551 

The spectrum of RX J09 11+0551 is presented in Fig. [T] and 
shows no definite detection of water maser emission within ± 2000 
kms -1 of the systemic velocity. The rms noise in the spectrum is 
1.6 mjy for a 10 kms -1 channel. This corresponds to a luminosity 
limit of < 7200 Lq (3 cr) when corrected for the estimated magni- 
fication from the lens (~ 22). 

4.2 IRAS 10214+4724 

Our Effelsberg spectra of IRAS 10214+4724 are presented inFig.|2] 
and show two tentative detections of water maser emission. The 
first epoch of data that were taken with Effelsberg shows a clear 
5 a peak in emission at 77 ± 1 km s -1 relative to the systemic ve- 
locity of the quasar. The peak of the fitted Gaussian is 4.7 mjy and 
the Gaussian FWHM of the line is 13 ± 3 kms -1 . The integrated 
total intensity of the potential maser line is 0.066 Jy km s - , giv- 
ing an isotropic luminosity of 3100 Lq (for a lensing magnifica- 
tion of 50). The noise in the spectrum within a 10 kms -1 chan- 
nel is 1.0 mjy channel -1 . The follow-up observations of IRAS 
10214+4724 with Effelsberg on 2008 March 14 and 15 did not 
detect this emission feature, but did show another emission line 
at —81 ±2 kms -1 that was significant at 4 a. This line has a fit- 
ted Gaussian peak flux density of 4.9 mJy and a FWHM of 18 ± 3 
kms -1 . The integrated intensity of the line is 0.096 Jy kms -1 and 
the isotropic luminosity is 4500 Lq. However, this second emis- 
sion feature was not detected on both days during the second Ef- 
felsberg observing epoch. The rms noise of the spectra taken on 
March 15 and 17 are 1.2 and 0.6 mJy channel -1 , respectively, 
within 10 km s -1 wide channels. 



The EVLA continuum map of IRAS 10214+4724 at 6.7 GHz 
shows an emission peak of just 0.19 mjy beam - close to the ex- 
pected position of the lens system, which is at the 2.5 a level (see 
Fig.(3](. A low sub-mjy flux density for IRAS 10214+4724 at 6.7 
GHz is consis tent with the earlier observations of the system at cm - 
wavelengths teowan-Robinson et al.lll99ll : lLawrence et alii 19931) . 
The EVLA spectrum, extracted at the position of the three merging 
images is also presented in Fig. [2] We find a 2.5 a peak in emis- 
sion at 72 ± 1 km s -1 , which is similar to the velocity of the ten- 
tative emission line found from the first epoch of Effelsberg obser- 
vations. Unfortunately, the bandwidth of the EVLA data was not 
large enough to span the velocity of the other tentative detection 
made during the second epoch of Effelsberg observations. The can- 
didate water line found using the EVLA has a peak flux density of 
1.4 mJy, a Gaussian FWHM of 3 ± 1 km s -1 and an integrated to- 
tal intensity of 0.010 Jy kms -1 . This corresponds to an isotropic 
luminosity of just 450 Lq (for a lensing magnification of 50). 



Our spectra of IRAS 10214+4724 that were taken with Ef- 
felsberg and the EVLA provide at best only tentative detections 
of water vapour emission at redshift 2.285, even though the can- 
didate lines from Effelsberg are significant (~ 5 a peaks). This is 
because we have not been able to verify these tentative lines from 
multiple observations with Effelsberg, or independently with the 
EVLA. It is certainly possible that the emission seen in the Effels- 
berg spectra is genuine, but is variable over the 4 month period 
between the observations. Water maser lines have been found to 
vary in flux density and line width over the time-scale of months 
before. For example, in the powerful FRII quasar 3C403, several 
maser lines at di fferent velocities a re seen to strongly vary over 
months to years jTarchi et al]|2007h . What is harder to explain is 
the variation seen within a few days during the second Effelsberg 
observation. However, given that the maser emission will come 
from a compact region, scintillation w ithin our own galaxy (e.g. 
Vlemmings, Bignall & Diamond 2007) or microlensing within the 
lensing galaxy cannot be discounted. The EVLA spectrum shows 
an emission feature which is close in velocity to one of the ten- 
tative lines found with Effelsberg, but given the line's low signif- 
icance (2.5 a) and width (~ 1 channel), it is not conclusive that 
this is a real detection. Note that the feature in the EVLA spec- 
trum would be too weak to have been detected during the second 
epoch of Effelsberg observations. Therefore, we conclude that there 
is only a tentative detection of water vapour emission from IRAS 
10214+4724. 
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Figure 1. The Effelsberg spectra of RX J091 1+0551 (top), H 1413+116 (top-middle), MS 1512-cB58 (bottom-middle) and SMM J16359+6612 (bottom). 
The spectra have a resolution of 10 km s -1 channel - 1 with an rms of 1.6, 0.9, 1.2 and 0.9 mJy channel" 1 , respectively for the four target lens systems. All 
velocities are relative to the systemic velocities of the quasars, and are in the heliocentric frame. The uncertainty in the optical systemic velocity of each target 
is around ± 100 kms -1 . These targets were observed using 100 MHz bandwidth. 



4.3 H 1413+117 - The Cloverleaf 



The spectrum of the H 1413+1 17 lens system is presented in Fig.Q] 
and shows no evidence of water maser emission. There was sig- 
nificant local RFI detected with the Effelsberg radio telescope be- 
tween 6.219 and 6.231 GHz, which resulted in part of the spec- 
trum between 600 and 1200 kms -1 from the systemic velocity of 
H 1413+117 having to be removed. The rms noise of the spec- 
trum within a 10 kms -1 wide channel is 0.9 mJy. This corre- 
sponds to an upper limit of < 7700 Lq (3 a) for an unlensed wa- 
ter maser (assuming a lens magnification of ~ 10). We note that 
our observat i ons ar e about half as sensitive as those presented by 
IWilneretal.Nl999h , but probe a much larger velocity range, out 
to around ±2000 kms -1 . The region of the Effelsberg spectrum 
affected by RFI was also not investigated by IWilner et alj dl999h 
because their spectrum was confined to velocities between —633 
and +127 km s ~~ 1 , relative to our definition of the systemic velocity 
of the quasar. Therefore, we also find no evidence of time-variable 
emission from this system. 



4.4 MS 1512-cB58 

The spectrum of MS 1512~cB58 is presented in Fig.Q]and shows 
no evidence of clear water maser emission within ± 2000 kms - 
from the systemic velocity. The rms noise of the spectrum is 
1.2 mJy for a 10 km s _1 wide channel. This corresponds to an up- 
per limit of < 2300 Lq (3 a) for the luminosity of any water maser 
emission within a 10 kms -1 wide channel. This calculation as- 
sumes a lensing magnification of 50. 



4.5 SMM J16359+6612 

The spectrum of our final candidate water maser galaxy is also 
shown in Fig. [T] We find no water maser emission in the spec- 
trum, which has an rms noise of 0.9mJy within a lOkms -1 wide 
channel. The upper-limit to the luminosity of water maser emis- 
sion is therefore < 1200 Lq (3 a) after correcting for the magni- 
fication from the cluster lens (~45). Comparing with the recent 
observations of SMM J16 359+6612 bv lEdmonds et alj d2009T) and 
Wagg & Momjian] d2009h . we find that our spectrum is about half 
as sensitive as those taken with the EVLA. The Effelsberg obser- 
vation confirms that there is no time-variable emission from this 
system and probes a much larger range in velocity. The EVLA ob- 
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Figure 2. The Effelsberg and EVLA spectra of IRAS 10214+4724. Top 
panel: the Effelsberg spectrum taken on 2007 November 12 shows a clear 
peak at +77 kms -1 from systemic. The rms is 1.0 mjy channel -1 (10 
kms -1 channel width). Middle panel: the Effelsberg spectra taken on 2008 
March 14 and 17. There is evidence of apeak at —81 kms -1 on March 14, 
but nothing was found on March 17. The rms is 1.2 and 0.6 mJy channel -1 , 
respectively. Bottom panel: the EVLA spectrum taken on 2008 May 30. 
Again, a faint emission line is detected at ~ 72 kms -1 . The rms is 
0.6 mJy channel -1 (4.3 kms -1 channel width). All velocities are rela- 
tive to the systemic velocity of IRAS 10214+4724 at redshift 2.285, and 
are in the heliocentric frame. The bandwidth of the Effelsberg and EVLA 
observations was 20 and 6.25 MHz, respectively. 



servations were processed with the old VLA correlator and gave 



information between —406 and 
nition of the systemic velocity. 



106 km s , relative to our defi- 



5 DISCUSSION 

In this section, we discuss our results and present an analysis of 
the water maser luminosity function that is based on our survey of 
gravitationally lensed quasars and star-forming galaxies. We also 
discuss the prospects for finding non-lensed water maser galaxies 
in the distant Universe with the next generation of radio telescopes. 



5.1 Water masers at high redshift 

There have been several attempts to detect water vapour from grav- 
itationally lensed galaxies. So far (including this work), eight grav- 
itationally lensed quasars and star-forming galaxies at high red- 
shift have been observed for 22.2 GHz water maser emission, with 
only one confirm ed detection (MG J0414+0534 at redshift 2.64; 
Impellizzeri et al. 2008). We present in Table [2] the eight systems 



40" 



30" 



20" 



10" 



Q 47°09'00' 



50" 




4C" 



10 24 m 37 s 36 s 35 s 34 s 33 s 
Right Ascension (J2000) 



32 a 



Figure 3. The continuum image of IRAS 10214+4724 taken with the EVLA 
at 6.7 GHz. An emission peak (0.2 mJy beam - 1 ; 2.5 a) is found ~ 4 arcsec 
from the position of the lens system (shown by the cross at 10 h 24 m 34.56 s 
+47°09'09.7"). Contours are shown at (-2, 2) X 0.08 mjy, the rms map 
noise. The restoring beam is shown as the grey ellipse in the lower left- 
hand side of the map. 



that have been observed and the 3 cr luminosity limits for those sys- 
tems where no water masers were found. With the exception of 
APM 08279+5255 at redshift 3.911, all of the lens systems have 
similar redshifts, forming a sample of galaxies at redshift ~ 2.5. 
Also, we find that for this sample, all of those systems with non- 
detections have luminosity limits which are below the measured lu- 
minosity of the water maser from MG J0414+0534. Therefore, we 
find the detection rate for powerful water masers with luminosities 
of ~ 10 4 Lq from galaxies at redshift ~ 2.5 to be one in seven. 

We focussed our sample selection on those galaxies with ev- 
idence of large reservoirs of dust through previous detections of 
CO and/or from a high FIR-luminosity to investigate whether dusty 
galaxies were more likely to host powerful water maser systems. 
Note that the estimated cold-dust masses of the sample are of or- 
der ~ 10 8 M fl and that MG 10414+0 534 has a cold-dust mass of 
~ 10 7 M Q |b arvainis & Ivisonl l2002). Our non-detections would 
suggest that the presence of dust is not enough to increase the like- 
lihood of finding water masers at high re dshift, in agreement with 
the findings of Wagg & Momiiarj ^20091) . A correlation between 
dust and water maser emission has b een tentatively found in the 
local Universe JCastangia etai1l2008h ; although, there is only one 
known local ultra -luminous infrared galaxy with water maser emis- 
sion (UGC5101: [Zhang et alJl2006T) . This correlation is strongest 
for water maser galaxies in the kilo-maser regime (^ 10 Lq), that 
is, those which have some or m ost of their water maser emission 
associated with star formation dSurcis et alj|2009h . Thus far, there 
is no evidence of such a correlation at moderate or high redshifts 
where galaxies are expected to have significantly more dust. CO 
may also not be a particularly good tracer of water maser emis- 
sion. Although CO does trace the cold gas of a galaxy, and in 
the case of MG J0414+0534 the water maser emission and the 
CO have consistent velocities, CO tends to trace lower density 
gas (10 4 to 10 6 cm -3 ) than is required for water maser emission 



(10 7 to 10 



11 cm 3 ). This implies that water maser emission and 
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the overall dust and gas masses of a galaxy are not related in a 
simple way. This then leads to the possibility that jets may be a 
key ingredient for luminous maser emission at high redshift. Rel- 
ativistic jets would provide a mechanism to chemically enrich the 
gas with water molecules through shocks, to collisionally excite 
the water molecules, as well as giving a strong continuum source 
to stimulate coherent emission. The next stage of our investiga- 
tion is to carry out a larger survey of high-redshift AGN that are 
gravitationally lensed, in particular, focussing on those systems 
with known radio jets. From surveys at radio wavelengths, there 
are around 3 5 gravitationally lensed AGN with known relativis- 
tic jets (e.g. Walsh, Carswell, & Wevmann 1979; Lawrence et al. 
ll984IWinnetai . 2000; Browne et al. 2003; Myers et al. 20"03b. 



5.2 Evolution of the water maser luminosity function 

The water maser luminosity function describes the co-moving num- 
ber density of water maser galaxies with luminosity, Lh 2 o, per log- 
arithmic interval in Lh 2 o and can be used to investigate the density 
and luminosity evolution of water maser galaxies as a function of 
redshift. The water maser luminosity function for the local Uni- 
verse was c alculated using up to 78 galaxies at distances less than 
235 Mpc bv lHenkel et alM2005h andlBennert et al.l d2009h . and was 
found to be oc L~^ ±0 ^. iBennert et alj d2009h also consid- 
ered the detection of SDSS J080430.99+360718.1 at redshift 0.66, 
and found that, given the sensitivity of their survey of 274 Seyfert 2 
galaxies with the Green Bank and Effelsberg telescopes, the detec- 
tion was consistent with the local water maser luminosity function. 

We have investigated whether the detection of a water maser 
with an unlensed luminosity of 8300 Lq from MG J0414+0534 at 
redshift 2.64 is consistent with the water maser luminosity func- 
tion f ound at low and moderate redshifts (see llmpelhzzeri et al.l 
2008). The probability of detecting the water maser from MG 
J0 4 14+0534 was estim ated by integrating the luminosity function 
of iHenkel et alj J2005h . with a slope of —1.4, over all luminosities 
greater than 10 Lq. It was assumed that only those water masers 
with luminosities stronger than 10 Lq are associated with AGN, 
which is the case in the local Universe. However, this split in lu- 
minosity may not be so clear at higher redshifts where it could be 
harder to separate the water maser emission associated with either 
star-formation or AGN activity. We found that for galaxies with 
water masers, the fraction that have maser emission greater than 
8300 I/© is predicted to be just 8 x 10 . This would imply that to 
find such a water maser system from a single pointing is rare, and 
would require a large volume to be surveyed for a detection to be 
made. However, we have carried out a targeted survey of AGN and 
star-forming galaxies. In the local Universe, the dete ction rate of 
water maser galaxies is approximately a few per cent teraatz et al.1 
2008). The combined probability of observing an AGN and find- 
ing an 8300 Lq water maser galaxy is estimated to be just ~10 . 
Therefore, our detection of at least one water maser galaxy at red- 
shift ~ 2.5 implies that there has to be some change in the water 
maser luminosity function from what is found out to moderate red- 
shifts (z ~ 0.7). The natural consequence of this evolution is that 
powerful water maser galaxies must be significantly more abundant 
at redshift ~ 2.5 when compared to the number found locally. As 
we have only a single detection at high redshift it is not possible to 
place any strong constraints on what form this change in the lumi- 
nosity function could be. For example, there could be a change in 
the slope of the luminosity function, or strong density and/or lumi- 
nosity evolution. 



5.3 Predictions for the SKA and pathfinders 

The next generation of radio telescopes, such as the EVLA, 
MeerKAT (Karoo Array Telescope) and the SKA (Square Kilome- 
tre Array), will allow deep and/or wide-field surveys to be carried 
out at frequencies where water maser galaxies can be found at cos- 
mologically interesting epochs (z > 1). To investigate the feasibil- 
ity of blind surveys with these telescopes, we have calculated the 
number of water maser galaxies we would expect to find. We give 
the expected specifications of the EVLA, MeerKAT and the SKA 
in Table[3]that were used for our calculations. 

We predicted the number of water maser gala xies by mig- 
rating the local water maser luminosity function of Henk el et al.1 
20051) . 



$ = (1 x 10 



-3 



L^ 2 o Mpc -3 (0.5 dex) -1 



(2) 



(where Lh 2 o is expressed in solar luminosities) between 10 Lq ^ 
Lh 2 o ^ 10 4 ' 66 Lq (i.e. twice the most luminous water maser 
known) and over the volume of sky that can be surveyed (defined 
by the field-of-view, frequency coverage and total instantaneous 
bandwidth). We assume that the luminosity function is a power-law 
without a knee before our luminosity cut-off. We also parameterize 
any evolution of the luminosity function with redshift using, 



(1 + *) 



(3) 



where m is 0, 4, and 8 for no, moderate and extreme evolution, 
respectively, up to redshift 2.2, after which there is assumed to be 
constant density evolution. 

How far we can probe down the luminosity function at differ- 
ent redshifts is dependent on our detection limit. This is calculated 
by first determining the system equivalent flux density (SEFD) us- 
ing, 



SEFD = 2 k T sys /A eS 



(4) 



where k is the Boltzmann constant and A c ff/T sys is given in Ta- 
ble[3] The limiting flux-density for the detection of a maser line is 
calculated using, 



SEFD 



(5) 



where r\ c is the correlator efficiency and is taken as 0.7, n po i is the 
number of polarizations, Av is the line width and t is the observing 
time. For our calculations, we assume a line width of 10 km s _1 and 
dual polarization observations. The luminosity limit is then found 
from Equation [T] 

We show the results in Table [4] of our simulations for a single 
pointing, 100 pointings and quarter-sky surveys (where appropri- 
ate) for the EVLA, MeerKAT and the SKA. We considered total 
integration times of 3000 h on-source to complete the surveys. The 
time per pointing is simply the total time divided by the number 
of pointings and frequency set-ups needed to cover the full range 
in frequency. We show the number of detections that are expected 
for the three evolutionary scenarios in Table [4] but here we only 
discuss the results for moderate density evolution, (1 + z) 4 . 

Taking the EVLA and MeerKAT results first, we find that 
these two instruments are not particularly useful for blind surveys 
of high redshift water maser galaxies, given their limited sensi- 
tivity and fields-of-view. However, the sensitivities of these tele- 
scopes are good enough to detect water masers out to redshift 1- 
2 with ~ 1000-5500 Lq using 15-30 h on-source (i.e. unlensed 
water masers with luminosities similar to or lower than that of 
MG J0414+0534). Therefore, the EVLA and MeerKAT will be im- 
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Table 2. Summary of gravitationally lensed quasars and star-forming galaxies that have been searched for 22.2 GHz water masers. All luminosity limits for 
non-detections have been calculated at the 3 a level and assume a 10 kms - 1 channel. Given the uncertainty in the magnification of any potential water maser 
systems, we also quote the magnification used for our calculations. 



Lens system 


Redshift 


Magnification 


Velocity coverage 


Luminosity 


Reference 








(kms^ 1 ) 


(£©) 




JV1G J0414+0534 


2.639 


35 


780 


< 3300 


Wiiner et al. ( 1999) 








i son 

1 JUU 


S'UVl 
ojUU 


iniijeiiizzen ei di. \z\j\Jo) 


APM 08279+5255 


3.911 


7 


760 


< 20000 


Ivison (2006) 


RXJ091 1+0551 


2.793 


22 


4200 


<7200 


this paper 


IRAS 10214+4724 


2.285 


50 


500 


<980 


this paper 


SMMJ1401 1+0252 


2.565 


4 


300 


<2900 


Wage & Momiian (2009) 


H1413+117 


2.561 


10 


780 


<4600 


Wiiner et al. (1999) 








4000 


<7700 


this paper 


MS1512-cB58 


2.720 


50 


4200 


<2300 


this paper 


SMM J16359+6612 


2.516 


45 


600 


< 1800 


Edmonds et al. (2009) 








300 


<650 


Waea & Momiian (2009) 








3600 


< 1700 


this paper 



portant for targeted searches of candidate water maser galaxies in 
the future, similar to those carried out here for the gravitationally 
lensed water maser project. We find that the proposed sensitivity 
and field-of-view of the SKA Phase 2 will be sufficient to detect 
many water maser galaxies out to high redshifts. In Fig.|4]the num- 
ber of predicted detections as a function of redshift for a single 
pointing, 100 pointings and quarter-sky surveys with the SKA are 
shown. The quarter-sky survey is predicted to detect at least 23500 
water maser galaxies out to around redshift 3, beyond which the 
sensitivity is not good enough to make detections in meaningful 
numbers. The time on-source per pointing for this wide-field sur- 
vey is around 5 s. The 100 and single pointing surveys are expected 
to have the sensitivity to detect around 3700 and 900 water maser 
galaxies, respectively, out to redshift 4.56, the redshift cut-off of 
our proposed surveys. 

These predictions are based on a number of assumptions, 
which we now discuss. First, we assume that the luminosity func- 
tio n of water ma s er gal axies is consistent with th at derived locally 
by iHenkel et all d2005l) and iBennert et al] d2009l) . This luminos- 
ity function is based on only 78 nearby galaxies, and does not 
probe the high luminosity regime (> 10 4 L ). We have adopted 
a high luminosity cut-off at 10 4 ' 66 Lq for making our predictions, 
but the resulting number of detections is not particularly sensitive 
to the choice of th is value. There is evidence for a knee in the 
Henk eTet al.l d2005T) luminosity function around 10 3 L@. However, 
this is thought to be due to there being few high luminosity masers 
in the limited volume surveyed, as opposed to a genuine steepen- 
ing of the spatial distribution of high luminosity water masers. We 
have tried to limit the effect of a potential cut-off in the luminosity 
function by only integrating the luminosity function up to twice the 
luminosity of the most luminous water maser known. Also, from 
the detection of a luminous water maser system at redshift 2.64, 
it was found that there must be some change in the local lumi- 
nosity function with redshift. Therefore, our predictions must give 
a lower limit to the number of water masers that we would ex- 
pect to find. The second major assumption in our predictions is that 
the number of potential host galaxies with the physical conditions 
needed for water maser emission will increase with (1 + z) 4 , that 
is, similar to the merger-rate of galaxies found from optical surveys 
dLe Fevre et ai]|200ol ; iRvan et alj 2008). We have tried to account 
for this uncertainty by using no and extreme evolution in the lumi- 
nosity function. Although no change in the number density of wa- 
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Redshift 

Figure 4. The number of water maser galaxies that are expected to be found 
(3 a level) with the SKA Phase 2. In total, 3000 h is used for either a single 
pointing, 100 pointings or quarter sky surveys. These calculations assume 
that the form of the luminosity function remains unchanged from that de- 
rived locally and the population of galaxies that host water masers evolve 
with (1 + z) 4 , up to redshift 2.2., after which there is constant density evo- 
lution. The error bars are at ± 2 a confidence level. 



ter maser galaxies with redshift is unlikely, extreme evolution may 
be possible. For example, in the case of optically b right quasars 
the evolution goes as (1 4- z) e up to redshift ~ 2.5 dBriggs|[l998l : 
iHewett, Foltz & Chaffedl 19931 ; ISchmidt et alll995h . 



6 SUMMARY 

We have presented new radio spectra for five gravitationally lensed 
quasars and star-forming galaxies at redshifts between 2.29 and 
2.80. We searched these galaxies for 22 GHz (rest frame) water 
masers that potentially could be used to study the properties of their 
central supermassive black hole. No new water masers were found, 
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Table 3. The specifications for the different arrays used to calculate the expected number of water maser galaxies from z = 0.5 up to z = 4.5 that can be 
found from blind surveys. 



Telescope 


Frequency 
(GHz) 


Number of 
antennas 


Antenna diameter 
(m) 


Aperture efficiency 


-lays 
(K) 


Total simultaneous bandwidth 
(GHz) 


(m 2 / K) 


EVLA 


4-12 


27 


25 


0.70 


38 


2.0 


245 


MeerKAT Phase 4 


10-14 


80 


12 


0.70 


30 


4.0 


210 


SKA Phase 2 


4-10 


2000 


15 


0.70 


30 


2.0 


8250 



Table 4. The predicted number of water maser galaxies found from blind surveys. The sky-areas are quoted at 10 GHz, except for the EVLA observations 
between 4 and 8 GHz, which are defined at 5 GHz. The luminosity limits are for 10 kms" 1 channels within the range of redshifts that are observable. The 
number of detections are for different evolutionary models given by (1 + z) m , where m is 0, 4, and 8 for no, moderate and extreme evolution, respectively. 
The L3 u denotes the 3 a limiting luminosity that can be detected across the full frequency range for the given time on-source per pointing. 



Telescope 


Frequency 
(GHz) 


Redshift 


Total time 
(h) 


Pointings 


Number of 
frequency settings 


Sky area 
(sr) 




(Lq) 


m = 


Detections 
m = 4 


m = 8 


EVLA 


8-12 


0.85-1.78 


3000 


1 


2 


1.7 x 10" 


-6 


90^150 








30 








3000 


100 


2 


1.7 x 10- 


-4 


900-4500 








100 




4-8 


1.78^1.56 


3000 


1 


2 


6.7 x 10~ 


-6 


450-3100 








250 








3000 


100 


2 


6.7 x 10" 


-4 


4500-31000 








850 


MeerKAT Phase 4 


10-14 


0.59-1.22 


3000 


1 


1 


7.3 x 10- 


-6 


30-160 








30 








3000 


100 


1 


7.3 x 10- 


-4 


300-1600 





10 


130 


SKA Phase 2 


4-10 


1.22^1.56 


3000 


1 


3 


4.7 x 10- 


-6 


7-110 


10 


900 


82000 








3000 


100 


3 


4.7 X 10- 


-4 


70-1100 


50 


3700 


330000 








3000 


6.7 X 10 5 


3 


3.14 




5600-90000 


400 


23500 


1800000 



although, we do have a tentative detection of water vapour from the 
quasar IRAS 10214+4724 at redshift 2.285. Our survey has now 
searched for water masers in six gravitationally lensed quasars and 
star-forming galaxies, doubling the number of high redshift lensed 
galaxies surveyed so far, and finding the onl y known high redshift 
water maser system (impellizzeri et al. 200j|). 

We have presented an analysis of the high redshift water maser 
lumi nosity function tha t is b ased on the results from this paper and 
from llmpellizzeri et alJd2008h . We again confirm that there must be 
some evolution in the water maser luminosity function to explain 
the detection of the luminous (8300 Lq) water maser system from 
MG J0414+0534 at redshift 2.64. We therefore conclude that wa- 
ter masers are more abundant at higher redshifts than found locally, 
and surveys with the next generation of radio arrays have the poten- 
tial to find many more water maser systems in the future. However, 
blind surveys for water maser galaxies will only be possible with a 
sensitivity like that of the proposed SKA (Phase 2). 
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